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^" ■ 1 INTRODUCTION 

°. 

Recent advances in our understanding of critical phenomena due to the application 
of the renormalization group (RG) approach 2 are now well known (see e.g. the text- 
books 3 ' 4 ' 5 ). The scale invariance at the critical point and the universality of certain 
features of critical phenomena can be explained by an RG transformation, and lead to 
a theory which provides a quantitative description of the critical behaviour of various 
£h ! thermodynamic quantities of interest. 

Specifically, the critical point in the RG " language" (in the context of our lecture 
we specify it as an equilibrium second order phase transition point) corresponds to 
the stable fixed point of the RG transformation, where the system is scale invariant. 
Asymptotic properties of the system are governed solely by the coordinate of the stable 
fixed point whereas non-asymptotic ones are defined in the region of the approach to 
the fixed point. In addition, the RG approach may permit the order of the phase 
transition occurring in a certain model to be determined. This is explored by studying 
the stability of the fixed points of the RG transformation: an absence of a stable fixed 
point is interpreted as an evidence of a fluctuation-induced first-order phase transition 
in the system under consideration. For models with no exact solutions or rigorous proofs 
of the existence of a second order phase transition, that is to say, for the majority of 
realistic models in statistical physics, RG provides a tool to check the order of the 
transition. 

Now let us define the kind of problem we are going to discuss in this lecture. It 
can be formulated as: What is the order of the normal-to-superconducting phase tran- 
sition?. According to the Bardeen-Cooper-Schriefer theory of the superconductivity, 
the normal-to-superconducting (NS) phase transition is a classical second order phase 
transition described by the Landau- Ginsburg Hamiltonian with a complex order pa- 
rameter corresponding to the wave function of the Cooper pairs. Taking into account 
the fluctuations of the order parameter one can find values of corresponding critical ex- 
ponents which in this case will coincide with the critical exponents of 0(n) symmetric 
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field-theoretical model for the case n = 2 (the XY model). Consequently, this leads 
to the assertion that the NS phase transition is described by the same set of critical 
exponents as the phase transition in the normal-to-superfTuid liquid. The latter set 
of exponents has been measured to a high accuracy 6 and also calculated by various 
methods 7 ' 8 < 9 > 10 < n < 12 < 13 > 14 . 

Taking into account that the corresponding "superfluid liquid" is charged in case 
of the NS transition, however, complicates the problem. This has been considered first 
of all by B. I. Halperin, T. C. Lubensky and S. Ma 15 and since then different ways of 
tackling it have been suggested. We shall discuss some of them briefly in the subsequent 
section. 

From an experimental viewpoint, when the question of the order of the supercon- 
ducting phase transition has been first discussed it has been considered more or less 
academic since, due to the large correlation length £ ~ 10 3 A, the first order character- 
istics can be seen only very near the phase transition; otherwise mean field behaviour 
is to be expected. This situation has changed after the discovery of high-T c supercon- 
ductors with correlation lengths of the order of lattice distances (£o ~ 1 A) 16 . Since 
then, critical effects have been observed in several experiments 11 ' 18 ' 19 ' 20 ' 21 ' 22, 23 ' 24 . 

Our main result presented in this lecture is that within the framework of the RG 
method applied to the original superconductor model, minimally coupled to the gauge 
field 15 , one still can demonstrate the existence of a second order phase transition with 
critical exponents distinct from those of a superfluid liquid. To prove this we shall 
consider the two-loop renormalization group functions for the model, paying particular 
attention to the fact that the loop expansion is asymptotic 27 ' 28 ' 29 . In this way we find 
several fixed points with new scaling exponents and a rich crossover behavior. Some of 
our results have been previously published in 30 ' 31 ' 32 . 

The lecture is, therefore, organized as follows. In the following section we shall 
give a brief review of the methods used to study the problem we are interested in. Then 
we describe the model of a superconductor, provide results of its study by the mean 
field approximation, obtain the expressions for the renormalization group functions in 
a two loop approximation and describe the results obtained subsequently without ap- 
plying any resummation procedure. Then in the next section we shall discuss several 
well-established examples in the modern theory of critical phenomena where the re- 
summation of the asymptotic series is applied. After that we shall present the pivotal 
section of our investigation: it is devoted to a study of the RG functions and of the cor- 
responding flows on the basis of the resummation technique. We obtain as a result that 
a stable fixed point is present, which is an evidence of a second order phase transition 
in the model. In the remainder of the paper the asymptotic and effective values for the 
critical exponents are calculated and we give expressions for the amplitude ratios. To 
conclude, these results are discussed in the closing section. 



2 NORMAL-TO-SUPERCONDUCTING TRANSITION: 1st OR 2nd OR- 
DER? 

The order of a phase transition may have severe effects on the physical quantities 
of a material. This is illustrated by the first-order liquid-gas transition phenomena 
like overheating and undercooling connected with the metastability at the transition. 
For second order phase transitions, divergences in physical quantities occur (in the 
thermodynamic limit, of course) leading to a dramatic increase of the specific heat or the 
scattering of the light (the critical opalescence) near the liquid-gas critical point. Similar 
dramatic changes are associated with the phase transitions which have occurred in the 
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early stages of the universe where the questions discussed here for superconductors also 
find relevance (for a recent review see 33 ). 

As mentioned in the Introduction, the question of the NS phase transition order 
becomes complicated when one accounts for the fluctuations of the order parameter 
being coupled to the "gauge field" (the vector potential of the fluctuating magnetic 
field created by the Cooper pairs), the fluctuations of which also diverge at long dis- 
tances. Posed for the first time more than 20 years ago 15 , this problem has remained a 
challenging one in the physics of the superconductivity until now. 

The theoretical model of Halperin, Lubensky, and Ma describing the relevant crit- 
ical behaviour is the usual 0{n) symmetric 4 model with the n/2-component complex 
field coupled to a gauge field which describes the fluctuating magnetic field created 
by Cooper pairs. The answer obtained in 15 states that due to the coupling to the gauge 
field, in the mean field approximation a third order term appears in the free energy 
of the superconductor and the NS phase transition is of first order. This mean field 
analysis is appropriate for the type-I superconductors 34 where the fluctuations in the 
order parameter have no significant effect on the thermodynamics of the transition. 

The case of type-II superconductors is more complicated, however, because fluc- 
tuations cannot be neglected. Studying the problem within Wilson-Fisher recursion 
relations 35 in the first order of e it has been found 15 that a stable fixed point (neces- 
sary, but not sufficient for a second order phase transition) exists only for the order 
parameter components number n > 365.9, far exceeding the superconductor case n = 2. 
The crossover near the first order phase transition has been examined 26 and the ex- 
pression for the crossover function of the specific heat is given by the one loop order 
perturbation theory. 

The kinetics of fluctuations arising from vortex pairs in a superconductor has been 
studied by means of numerical simulations 36 . The result leads to the conclusion that 
there exist nucleation processes typical for the first order phase transition, confirming 
the mean field and RG results of 15 . Note, though, that the mean field analysis applied 
to the Ginzburg-Landau free energy of a superconductor 37 including a Chern-Simons 
term leads to quantitatively different behaviour: for different values of the topological 
mass, there occurs in a system either a fluctuation-induced first order phase transition 
or only a second-order transition. This result is also confirmed in the framework of 
one- loop RG calculations 37 . 

The occurrence of a first order phase transition has been also found in the massless 
scalar electrodynamics 38 ' 39 and confirmed to a linear order in e for the n-component 
Abelian Higgs models by an explicit construction of the coexistence curve and the 
equation of state 40 . In addition, the CP^ 1 non-linear sigma model, being related 
to the superconductor model in the limit of infinite charge by means of the (2 + e) 
expansion, has been shown 41 to posses a similar behaviour observed further by 15 . 

Results for type-II superconductors obtained by e -expansion methods 15, 26 ' 40 ap- 
pear to be stable against the influence of different physical factors such as the possibility 
of another (non-magnetic) ordering, the presence of a disorder and a crystal anisotropy 
when the study is performed by means of a strict expansion. The scaling behaviour 
of a superconducting system with an additional non-magnetic ordering studied by e- 
expansion methods provides one more example of a system where a weak first-order 
phase transition occurs 45 . The analysis of the influence of quenched impurities on the 
critical behaviour of superconductors when taking account of the magnetic field fluctu- 
ations demonstrates 46 ' 47 the appearance of a new stable fixed point for 1 < n < 366. 
It has been shown 47 , however, that this describes the critical behaviour in the range of 
space dimensionalities d c (n) < d < 4 with d c {2) = 3.8 and results in a first order phase 
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transition. 

The RG flow for the superconductor model with quenched impurities has been 
found 48 to exhibit a stable focus surrounded by an unstable limit cycle. The second 
order phase transition behaviour is found to show up inside the limit cycle. Introduc- 
ing random fields with short and long range correlations does not lead to a second 
order behaviour in the region of (d,n) near (3,2) 49 . Note however, that studies of 
the influence of quenched and annealed gauge fields on the spontaneous symmetry 
breaking, performed in terms of Helmholtz free energy 42 lead to the conclusion that in 
the first nontrivial or one-loop approximation in the annealed model the spontaneous 
symmetry breaking occurs through a first order transition for d = 2, 3 whereas the 
quenched model displays a continuous phase transition. A more complicated account 
of fluctuations in the annealed model changes the nature of the transition to a con- 
tinuous one, but the spontaneous symmetry breaking is absent in the model with a 
quenched disorder 42 . The combined influence of the crystal anisotropy, the magnetic 
fluctuations, and the quenched randomness on the critical behaviour of unconventional 
superconductors 50 studied by means of the RG analysis within the e-expansion 51 gives 
that only fluctuation-induced first order transitions should occur in unconventional 
superconductors in the vicinity of the critical point. 

The results of mean field calculations, however, have been questioned by Lovesey 42 , 
mentioned above, who has shown that taking into account the gauge field fluctuations 
in the calculation of the free energy leads back to a second order phase transition. 
A further indication of a second order phase transition has come several years later 
when this question has been studied as a lattice problem by means of Monte-Carlo 
calculations and duality arguments 43 . These results have confirmed that there are 
scenarios of the NS transition that differ from those obtained in 15 . Namely, the NS 
transition has been found to be of second order, asymptotically equivalent to that of a 
superfluid with a reversed temperature axis. Subsequent MC simulations 44 performed 
in different regions of couplings lead to the result that the NS transition is strongly 
first order deep in the type-I region and becomes more weakly first order moving in the 
direction of the type-II region. Beyond a certain point the data reported in 44 suggests 
a second-order transition. The corresponding 0{n) nonlinear a-model coupled to an 
Abelian gauge field studied near two dimensions by 2 + e expansion 53 does not show a 
first order phase transition either. 

The existence of a tricritical point, where the order of the phase transition changes 
from second into first, has been predicted 58 by representing the 3d superconductor 
model by a disordered field theory. The position of the tricritical point is located 
slightly in the type-I region for the values of the Ginzburg parameter 34 k < 0.8/\/2. 
Starting from the dual formulation of the Landau- Ginzburg theory, by means of the RG 
arguments, it has been shown that the critical exponents of the NS transition coincide 
with those of a superfluid transition with a reversed temperature axis 59 . But, while the 
correlation length critical exponent of the normal-to-superconducting transition is pre- 
dicted to coincide with the ordinary 3D XY model, the divergency of the renormalized 
penetration depth is characterised by the mean field value v = 1/2 59 . 

The influence of the critical fluctuations on the order of the NS transition has 
been reconsidered on the basis of the field theoretical RG ideas in 30 . Here the two- loop 
flow equations 30 for the static parameters and the (^-functions 52 are obtained and it 
is indicated that a stable fixed point possibly exists and, as a consequence, a second 
order phase transition may occur. An attractive feature of the flow found in 30 is that 
it discriminates between type-I and type-II superconductors, depending on the initial 
(background) values of the couplings. For small values of the ratio (coupling to the 
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gauge field) / (fourth order coupling) (appropriate for type-II superconductors) the flow 
comes very near to the fixed point of the uncharged model but ends in a new supercon- 
ducting fixed point. For large values of the ratio (type-I superconductors) the flow runs 
away. For intermediate values of the ratio, the critical behavior may be influenced by 
a second (unstable) superconducting fixed point with scaling exponents quite different 
from those for the uncharged model. 

A flow picture qualitatively similar to 30 has been obtained in 54 by solving the 
model of a charged superconductor approximately with the help of nonperturbative 
flow equations, a method which appears to give very encouraging results for critical 
scalar field theories 55 ' 56 . Depending on the relative strength of the ratio (coupling to 
the gauge field) / (fourth order coupling) a first or a second order phase transition has 
been found. An approximate description of the tricritical behaviour has been given as 
well as an estimate of the correlation length critical exponent v and the pair correlation 
function critical exponent 77, which give us the second order phase transition has been 
reported. Depending on three different assumptions for the stable fixed point value 
of the coupling to the gauge field, the following values are obtained in two successive 
truncations of the potential; (t],u) = [(-0.13, 0.50); (-0.20, 0.47)], [(-0.13, 0.53); (-0.17, 
0.58)], [(-0.13, 0.59); (-0.15, 0.62)], indicating that the critical exponents belong to the 
physical region i] > 2 — d and v > 0, independent of the truncation, clearly pointing 
towards a second order phase transition. 

In the context of baryogenesis the question of the NS phase transition order has 
been considered within the two loop approximation in 57 and the effective potential has 
been calculated. The e-expansion has been applied to the electroweak phase transition 
in order to estimate various parameters of it in leading and next-to-leading orders in 
e, including the scalar correlation length, the latent heat, the surface tension, the free 
energy difference, the bubble nucleation rate, and the baryon nonconservation rate. Of 
course, the transition has been found to be a first order since only run away flows occur 
in the strict e-expansion perturbation theory. Note, that in the electroweak scenario of 
the baryogenesis there exists a so-called Sakharov requirement which is met when the 
transition is strongly first order rather than second order. 

The NS transition problem has been also studied by an analytical method which 
is not based on e or 1/n expansions. Using a non-perturbative method of solving the 
approximate Dyson equation for arbitrary d and n 60 it has been found 61 that the NS 
phase transition is governed by a "charged" fixed point. The value of the pair correlation 
function critical exponent 77 at d = 3, n = 2 is r?(3, 2) = —0.38. It is interesting to note 
that although the result for rj appears to be a well-behaved function of d and n, it breaks 
down at the critical value n c ~ 18 when expanded in e. Hence, the conclusion is drawn 
that the results of the ^-expansion obtained in 15 and, in particular, the absence of a 
stable fixed point solution for n < n c ~ 365.9 are to be interpreted as the breakdown of 
the e-expansion rather than a fluctuation-induced first-order phase transition. On the 
other hand, near d = 4 the results of 61 are in a good agreement with the e-expansion 
data 15 for high n (n > 366). 

Recently, the same problem studied with the RG technique in a fixed dimension 
d = 3 within the one-loop approximation has given an evidence of an attractive charged 
fixed point distinct from that of a neutral superfluid, leading, in particular, to the 
correlation length critical exponents values v ~ 0.53 and rj ~ —0.70 62 . Considered 
in the form of the continuum dual theory 63 , however, the magnetic penetration depth 
has been shown to diverge with the XY exponent, contradicting results mentioned 
above 30 ' 54 ' 61 ' 62 . To investigate this discrepancy, MC simulations of the 3D isotropic 
lattice superconductor in a zero external magnetic field have been performed. This 
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results in the conclusion that there is a single diverging length scale consistent with 
the universality of the ordinary 3D XY model 65 . Further applications of the model 
containing coupling to the gauge field have been suggested in the context of the quantum 
Hall effect 66 . 

Now, let us give a brief resume of the experimental data relevant to our study. 
As mentioned in the Introduction, the effects of thermodynamic fluctuations are gen- 
erally small in conventional low-T c superconductors because of their low transition 
temperatures and large coherence length. In contrast, high transition temperatures 
and small coherence lengths mean that critical fluctuations are relevant in high-T c 
superconductors. Though critical fluctuations in high-T c superconductors have been 
observed in a series of experiments ( see e.g. 17 ' 18 ' 19 ' 20 ' 21 ' 22, 23 ' 24 ) their interpretation 
has been changed somewhat. Deviations from the mean field (i.e. first order) behaviour 
have accounted for either by 3d Gaussian fluctuations (giving, in particular, values for 
the specific heat critical exponent a and the correlation length critical exponent v. 
a = v = 1/2) 17 ' 19 or by a nontrivial XY behaviour characteristic for an uncharged 
superfluid (with v ~ 2/3 and logarithmic divergences in a) 18 ' 20 ' 22,24 . Measurements 
of the heat capacity 22, 24 , the magnetization and the electric conductivity 22 of single- 
crystal samples of YBa2Cu 3 7 ^ x in a magnetic fields near T c support the existence 
of a critical regime governed by the XY-like critical exponents 18 ' 20 ' 22 ' 24, 25 , a similar 
conclusion follows from the crossover analysis of the zero-field heat capacity on a com- 
parable sample 21 . The maximum applied magnetic field for which the 3D XY scaling 
is valid, however, differs for various materials 67 . 

To conclude this brief review it is worth mentioning one more physical interpre- 
tation of a charged field coupled to the gauge vector potential. Namely, this is the 
nematic-smectic-A transition in liquid crystals 68 ' 69, 70, 71 ' 72, 73 . The nematic phase 
is an orientationally ordered but translationally disordered phase, rodlike molecules 
are aligned with their long axes parallel to the director and the smectic-A phase con- 
tains layers of molecules with their long axes perpendicular to the layer. It has been 
proposed 68 ' 69 that this transition can be described by a model similar to those describ- 
ing the NS transition in the charged case 15 . Now the smectic order parameter (being a 
complex field \P(r) that specifies the amplitude and the phase of the density modulation 
induced by the layering) is coupled to the director fluctuations. Contrary to the NS 
transition, the nematic-smectic A transition is characterized by a critical region in the 
experimentally accessible range. For certain materials it has been indeed shown 72 that 
both the latent heat data obtained through an adiabatic scanning calorimetry as well 
as independent interface velocity measurements near the Landau tricritical point can 
be fitted by a crossover function consistent with a mean field free energy density that 
has a cubic term 15 , implying that the nematic-smectic-A transition is a weakly first 
order. Many liquid crystals, though, appear to exhibit a continuous nematic-smectic-A 
transition (see 73 and references therein). High- resolution heat-capacity and x-ray stud- 
ies of the nematic-smectic-A transition performed during the past twenty years (see 73 
for a comprehensive review) show complex systematic trends to crossover, from three- 
dimensional XY to tricritical behaviour and an anisotropic behaviour due to a coupling 
between the smectic order parameter and director fluctuations. 



3 THE MODEL AND ITS "NAIVE" ANALYSIS 

Now it is well-known, that the influence of the order parameter fluctuations on the 
NS transition can be described by the Landau- Ginsburg free energy functional: 
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m = J dM|l0o| 2 + ||(V0 O | 2 + ^|0o| 4 }, (1) 

to is temperature-dependent, uo is a coupling constant and the complex order parameter 
4>o is connected with the wave function of Cooper pairs. The Cooper pairs are charged 
and therefore create a fluctuating magnetic field which leads to the appearance of 
additional terms in the free energy functional. Note, that this is not the case of a 
normal-to-superfluid transition in a neutral (uncharged) fluid, which is well described by 
([l]) without any modification. Describing the fluctuating magnetic field B by the vector 
potential A (B = rot A) and adding to (]!]) the minimal coupling between the fluctuating 
vector potential and the order parameter one obtains the free energy functional F[&, A], 
originally considered in 15 for a generalized superconductor in d dimensions with a d- 
dimensional vector potential A and the order parameter \1/ consisting of n/2 complex 
components. 

One can now describe the fluctuation effects by an Abelian Higgs model with the 
gauge invariant Hamiltonian 15 : 



H = J d rf *{|l*o| 2 + ||(V- ze Ao)*oi 2 + ^>o| 4 + i(V x A ) 2 }, (2) 

which depends on the bare parameters to, eo, uq. The parameter to changes its sign at 
some temperature, the rest of the parameters are considered temperature- independent. 
When the coupling constant eo = no magnetic fluctuations are induced and the 
model reduces to the usual field theory (jl]) describing a second-order phase transition 
and corresponding in the particular case n = 2 to the superfluid transition in 4 He. 

The mean field results for the critical behaviour of the model with a free energy 
functional F[&, A] corresponding to the Hamiltonian @ have been already reported in 
the original paper of Halperin, Lubensky and Ma 15 . In the framework of the mean field 
theory one can determine that systems characterized by the free energy functionals F[<j)] 
([I]) and F[^f, A] (n = 2) possess a qualitatively different critical behaviour. Neglecting 
the order parameter fluctuations (in accordance with the Ginzburg criterion this may 
be done for a good type-I superconductor) shows that depending on the sign of t the 
free energy ([I]) is minimized by the value of the order parameter = for t > or by a 
non-zero value, when t < 0, and that the appearance of the non-zero order parameter 
is continious; in the system under consideration a second order phase transition occurs. 

When applied to the free energy functional F[ty, A], however, the mean field theory 
predicts a qualitatively different behaviour. Defining the effective free energy -F e fj[\l/] 
as a function of the single variable \1/ by taking the trace over the configurations of the 
vector potential A one finds 15 that the expression for F[^>] will contain a term which 
has a negative sign and is proportional to l^] 3 . Such a term inevitably leads to a first 
order transition; F e g [^f] develops a minimum at a finite value of \& when the coefficient 
of the quadratic term is still slightly positive. 

As it has been already mentioned, the above reasoning is appropriate for a type-I 
superconductor. The case of type-II superconductors is considerably more complicated. 
Here, fluctuations in \l/ cannot be neglected and one must choose an appropriate tech- 
nique to study the problem. Originally, the critical behaviour of the model (0) in the 
presence of order parameter fluctuations has been studied in 15 with the help of Wilson- 
Fisher recursion relations 35 in the first order of e = 4 — d and in result it is shown, 
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in particular, that the second order phase transition is absent for n = 2 in the re- 
gion of couplings appropriate for the type-II superconductor. We will reproduce these 
£-expansion results below and then continue to analyze the problem further. 

In order to describe the long-distance properties of the model ([|) arising in the 
vicinity of the phase transition point we shall use a field-theoretical RG approach. Two- 
loop results 30 for the RG functions corresponding to (|2|) are obtained on the basis of 
a dimensional regularization and a minimal subtraction scheme 74 , defining the renor- 
malized quantities so as to subtract all poles at £ = 4 — d = from the renormalized 
vertex functions. The renormalized fields, mass and couplings are introduced by: 

*o = 4 /2 ^, A = 4 /2 A, to - t 0c = Z t ZyH, 

el = Z\Z~^ e s 7/ 'Sj 1 , u = ZuZ^u^Sj 1 (3) 

with e — 4— d. Here, \x is an external momentum scale, to c is a shift which for the results 
considered here can be set to zero, and Sd stands for the surface of a ci-dimensional 
hypersphere: Sd = 2 l ~ d ix~ d l 2 /Y{d/2). The Z-factors are determined by the condition 
that all poles at e = are removed from the renormalized vertex functions. 

The RG equations are written bearing in mind the fact that the bare vertex func- 
tions r^' are calculated with the help of the bare Hamiltonian (0) as a sum of one- 
particle irreducible (1PI) diagrams 75 : 



r£' M ({r}, {R}) =< * (ri) • • • Mr N )MRi) ■ ■ ■ M^m) >ipi, (4) 

which do not depend on the scale /i and their derivatives with respect to fi at fixed 
bare parameters are equal to zero. So one gets 



9 r JV,Mi 9 ry-N/2 r7 -M/2- r N,Mi n / K \ 

'° = ^~d~ * A R '° = ^ ' 

where the index means a differentiation at fixed bare parameters. So, the RG equa- 
tions for the renormalized vertex function T^' M will be: 



where / = e 2 , = — (t an d the RG functions read 



Ou Of 
p u (u, f) = /j,—\ , (3 f (uJ) = fi—\ , 

dlnZv d\nZ A d\nZ t 

Using the method of characteristics the solution of the RG equation may be written 
formally as: 



T% M (t, u, /, //) = X(l) N / 2 (X'(l)) M ^ M (Y(l)t, u(l), f(l),»l), (8) 
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where the characteristics are the solutions of the ordinary differential equations: 



J^lnX(0 = C*(«(0>/(0)> lj l \nX\l) = UHl)J{l)), 

lj l \nY(l) = Uu(l)J(l)), (9) 

lj lU (l) = p u (u(l), /(Q), lj t f{l) = (3 f (u(l), /(/)) (fO) 



with 

x(f)=x'(f) = r(f) = i, u(i) = u, f(i) = f. (ii) 

For small values of I, the equation ([Bp is mapping the large length scales (the critical 
region) to the noncritical point 1 = 1. In this limit the scale-dependent values of the 
couplings u{l), f(l) will approach the stable fixed point, if it exists. 

The fixed points u*, f* of the differential equations ( |10|) are given by the solutions 
of the system of equations: 

tow n = o, 

p u ( u *, n = o. (12) 

The stable fixed point is defined as the fixed point where the stability matrix 

Ba = , Ui = {u, /} (13) 

possess positive eigenvalues (or if complex, the eigenvalues with positive real parts). 
The stable fixed point corresponds to the critical point of the system. As we have 
mentioned above, in the limit / — > (corresponding to the limit of an infinite correlation 
length) the renormalized couplings reach the values they have in the stable fixed point. 

Now we can write the results for the RG functions obtained in a two-loop ap- 
proximation 30 following the above described procedure in frames of the dimensional 
regularization and the minimal subtractions schemes. From a Ward identity one has 
= Z e , and the remaining Z-i actors are to be found from the corresponding vertex 
functions T 2 ' , T ' 2 , and T 4,0 . Since the gauge field is massless, the renormalization has 
been performed at a finite wave vector. The results in the two-loop order read: 

= 1 + i{3e 2 - u 2 {n + 2)/144 + e%n + 18)/4e - (lln + 18)/48]}, (14) 

Z K = 1 + ^{-ne 2 /6 - ne 4 /2}, (15) 

Z t = l + -{{n + 2)u/6 + u 2 [{n + 2)(n + 5)/36e - (n + 2)/24] + 

ue 2 [-(n + 2)(l/2e - 1/3)] + e 4 [(3n + 6)/2e + (5n + l)/4]}, (16) 

Z u = 1 + -{(n + 8)u/6 + 18e 4 /u + u 2 [{n + 8) 2 /3Ge - (5n + 22)/36] + 

ue 2 \-(n + 8)/2e + (n + 5)/3] + e 4 [(3n + 24)/e + (5n + 13)/2] + 

e 6 /«[3(n + 18)/e-7n/2-45]}. (17) 

Following the standard procedure one then shows that the expressions for (3- 
functions in the two-loop approximation will be: 
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71 

Bf = -£/ + g/ 2 + nf, (18) 

A = _ £U + ^_^±ii M 3_ 6M/ + 18/2 

+ ^„ V+ Z™„ /M7 „ + 90)/ 3. (19) 



The previous analysis of the equations of type (|I8"D , fll9|) either on a one-loop 15 or 
two-loop level 30 has been based on the direct solutions of the equation for the fixed 
point. In the present study we want to point out that the series have a zero radius of 
convergence and they are known to be asymptotic at best. Therefore, some additional 
mathematical methods have to be applied in order to obtain a reliable information on 
their basis. 

We shall start by recalling the results of an e 2 -expansion for /^-functions 15 ' 30 . For 
the second order in e one obtains three fixed points: Gaussian {u* G = f* G = 0), "Un- 
charged" (u* u ^ 0J* U = 0) and "Charged" (u* c ^ 0J* C ^ 0), to be denoted as G, 
U, C. The expressions for them read: 



G: u*° = 0, f* G = 0, (20) 
U: u* u = u^e + u^e 2 } f* u = 0, (21) 
C : u* c = v°e + u c 2 e\ f* c = fie + fie\ (22) 



where 



u_ 6 v _ 18(3n + 14) 



,c 



n + 8 A (n + 

3(n + 36) + (n 2 - 360n - 2160) 1/2 



XL 

1 3n{n + 8) 

a\ n \n J 

with 

n + 8 r 36 3n + 14/ r \ 3 r /6\ 3 

ai = l + — ai = ^-\u 1 J -enu^-J + 

36 /6 V (n + 5)4 / y _ 3(71n+174) + f6_V + gQ) _ 

Almost all physical results concerning the phase transition described by the field 
theory (0) are to some extent based on the information given by (^) - (p2|). The main 
ones read: 

(i) the fixed point U is unstable with respect to the presence of an /-symmetry at 
d < 4 with a stability exponent 



u 
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Figure 1. /3-functions of the model of superconductor u (u, /), 0f(u, /) in one-loop approximation 
for d = 3, n = 2. The fixed points have coordinates (u* = f* = 0), (u* = 0.6, /* = 0). 



(ii) the fixed point C appears to be complex for n < n c = 365.9 already 15 on a 
one-loop level. The stability exponent is given by 



and on the two-loop level it reads: 



A M — — es, s — 



1 + 



36\ 432(n + 8) 



n l/2 



n 



rr 



leading to an oscillatory flow of u in a one-loop order below n c with the soluti- 

26, 30. 



on 



f(l) 



6fV 



u(l) = f(l) 



6 + nef(l- 
n 



2(n + 8) 



<-iy 

s tan 



(23) 



lie 



In f(l)r L l 



/2(n + 8)u n + 36 

+ arctan H 

\ sn f ns t 

here / and u are the initial parameters for I — 1; 



n + 36 



n 



(24) 



(iii) from the condition of the positiveness of the fixed point coordinate /* (/ = e 2 ) it 
follows that at e — 1, n has to be larger that 36. This questions the applicability 
of the e-expansion for n = 2 to d = 3. 
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beta 




Figure 2. /^-functions of the model of superconductor f3 u (u, /), f3f{u, /) in two-loop approximation 
for d = 3, n = 2. Only the Gaussian fixed point u* = f* — survives. 



So, it follows that for the "superconductor" case n = 2 which is the most interesting 
from a physical point of view a stable fixed point does not exist and, therefore, the 
observed phase transition is of a first order. 

Now we will study the RG equations in the minimal subtraction scheme in the 
framework of d = 3 theory 11 ' 12 ' 13 putting e = 1 in expressions for the RG functions 
and studying the perturbation theory in powers of coupling constants. The latter 
corresponds to the number of loops in Feynman diagrams and one develops, therefore, 
the perturbation theory in a successive number of loops. Direct calculations based 
on equations (|i8|), ( |i9|) at fixed d = 3 do not bring qualitatively new features to the 
described above analysis. In the one-loop approximation, leaving square terms in (|18D , 
( |T§D one finds that only one nontrivial fixed point u* = 6/(n + 8), /* = exists. The 
/3-functions (3 u (u, f), Pf(u, f) in the one- loop approximation at d — 3, n — 2 are shown 
in the Fig. |T[ The simultaneous intersection of the surfaces corresponding to both 
functions with the plane (3 = results in the fixed points; for the n = 2 they have 
coordinates u* = f* = and u* = 0.6, /* = as seen in the picture. In the two-loop 
approximation only the Gaussian fixed point survives, as one may see from Fig. |^. 

Nevertheless one should note that such a straightforward interpretation of the 
above expansions data has been questioned and the way of analyzing the series for 
/3-functions (0), (0) avoiding the strict e-expansion and exploiting the information 
on the accurate solution for the pure model at d — 3, has been proposed 30 . Also, from 
the comparison of e-expansion data for /* (giving positive values of /* only for n > 36) 
with the value of /* obtained without the e-expansion (remaining positive for all n) it 
has been conjectured that the lower boundary for n resulting in the negative /* may 
be an artifact of the expansion procedure. Let us now consider expressions for the 
RG functions more carefully, paying attention to their possible asymptotic nature and 
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treating them with a resummation procedure. 



4 RESUMMATION 

The resummation technique appropriate for critical phenomena and applied to the 
asymptotic series for the RG functions enables one to obtain extremely accurate values 
of the critical exponents 76 . In fact, the asymptotic nature of the series for the RG func- 
tions has only been proven in the case of the 4 model containing one coupling of an 
0(n)-symmetry (the n- vector model). The high-order asymptotics for these series are 
known 27 ' 28 ' 29 in analytical form as well. These results give the possibility of obtaining 
precise values for the critical exponents of the n-vector model by a resummation of 
the corresponding series for the renormalization group functions (see e.g. 7 ' 8 ' 11 ). As 
far as we know, no information similar to that obtained in 27 ' 28 ' 29 for the "uncharged" 
case (/ = 0) is available for the "charged" model we are considering here. For mod- 
els containing several couplings of different symmetries the asymptotic nature of the 
corresponding series for the RG functions is a generally accepted belief rather than a 
proven fact. 

As it will be important in the course of our future analysis we shall mention here 
the weakly diluted n-vector model, describing the ferromagnetic ordering in a system 
of Ni classical n-component "spins" located in N sites of a lattice (Ni/N < 1) and 
quenched in a certain configuration. Using the replica trick in order to perform the 
quenched averaging one concludes 82 that an effective Hamiltonian of such a model 
contains two fourth order terms of a different symmetry and reads: 



h= fd^W[\^ a \ 2 +rn»\ 2 ]- v i Err +^E(ri 2 ) 2 }> (25) 

J Z a=l ' ° \a=l ) ^' a=l 

where <p a is an n-component vector <p a = ((j) 01 ' 1 , 4> a ' 2 , ■ ■ ■ , 4> a ' n ); Uq > 0,Vq > are 
bare coupling constants; m is the bare mass and in the final results a replica limit 
m — > is to be taken. The RG functions for this model are obtained in the form 
of a double series in renormalized couplings u, v and the asymptotic nature of the 
series has not been proven for this model up till now 84 . Nevertheless, the appropriate 
resummation technique (applied as if these series are asymptotic ones) enables one to 
obtain accurate values for critical exponents in three dimensions 77 ' 78 ' 79 ' 80 ' 81 and to 
describe (in the n = 1 case) the experimentally observed crossover to a new type of a 
critical behavior caused by the weak dilution 85, 86 . These results have been confirmed 
by Monte-Carlo 87 ' 88 and Monte-Carlo RG 89 calculations. 

The two main ways of a resummation commonly used for the asymptotic series 
arising in the RG approach are: (i) a resummation based on the conformal mapping 
technique and (ii) the Pade-Borel resummation. Case (i) is based on the conformal 
transformation, which maps part of the analytical domain containing the real positive 
axis onto a circle centered at the origin and the asymptotic expansion for a certain 
function is thus re- written in the form of a new series (see 8 ). This resummation, 
however, is based on the knowledge of subtle details of asymptotics (the location of 
the pole, the high-order behavior) which are not available in our case. 

In the absence of knowledge about the singularities of the series, the most appro- 
priate method which can be used to perform the analytical continuation is the Pade 
approximation resulting in the Pade-Borel resummation technique (ii) (see e.g. 7 ). We 
are going to apply it to the special case of / = 0, so let us disucss it in detail. 
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Starting from the Taylor series for the function f(u): 



f(u) = J2 c i uj > ( 26 ) 

j>0 

one constructs the Borel-Leroy transform 

^gr^TT)^' (27) 

with r(x) the Euler's gamma-function and p some positive number 90 . Then one repre- 
sents ( PTD in the form of the Pade approximant ^^(wt): 

L 

Ff L a U x ) = 1 ir— (28) 

(in the subsequent analysis, proceeding in a two-loop approximation we will use the 
[1/1] Pade approximant) and the resummed function will be given by: 

oo 

f Res (u) = Jdt e - * t p F Ps * d6 (ut). (29) 
o 

The resummation scheme fl2"T|) - ( p9|) of the (asymptotic) series in one variable 



(|26|) is easily generalized to the two-variable case when the series is given in a form: 



f(u,v)= Ciju'v*, (30) 
jJ>o 



with the Borel-Leroy transform 



nu,v,t ) =z WT f^ TT) (3D 

The procedure is aimed to help how to choose an appropriate form of the analytic 
continuation of the series fl3"T|) . Two most common methods to proceed are the Borel 



resummation combined with Chisholm approximants and the Borel resummation of the 
resolvent series, presented in a form of the Pade approximant. For the first method, in 
order to write an analytic continuation of the series ([Jl]) one uses the rational approx- 
imants of two variables: so-called Canterbury approximants or a generalized Chisholm 
approximants 91 ' 93 which are a generalization of Pade approximants in the case of sev- 
eral variables, representing ( [31] ) in a form: 

^Chisholm/ ^ = ij_ ( 32 ) 

M 

(sums in the numerator and denominator are limited by the condition of the cor- 
respondence between known numbers of terms in the initial series and that in the 
approximant). Again, the resummed function is given by an integral (f29|): 

oo 

/ Res (u, v)= Jdt e _< t p F Chisholm (ut). (33) 
o 
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Figure 3. /3-functions of the diluted Ising model f3 u (u,v), (3 v (u,v) in one-loop approximation for 
d = 3. The fixed points have coordinates (u* = v* = 0), (it* = 0.667, u* =0). 



Proceeding with the second method, one writes for the series of two variables 
the so-called resolvent series T{u, v, r) 92 ' 93 introducing an auxiliary variable r, which 
allows the separation of contributions from different orders of the perturbation theory 
in the variables u, v: 

F{u,v,t) = Y, CijiurYivry, (34) 

i,j>0 

f(u,v) = F(U,V,T= 1). 

Now the resummation of the series J-(u, v, r) is performed with respect to the variable 
t as for the series in a single variable, applying the above described scheme ( p7|) - (p9|). 



Let us illustrate how the resummation procedure works in the case of the effective 
Hamiltonian fl25|). In order to make a direct comparison with the superconductor case, 



let us take the /3-functions obtained for the model fl25|) in the minimal subtraction 
scheme in the two-loop approximation, though high-order results are available for this 
model 81, 95 as well as results 77 ' 78 ' 79, 80 obtained in the d = 3 massive field theoretical 
approach 94 . The expressions for the /3-functions, corresponding to the renormalized 
couplings u, v in the replica limit m — > for the Ising model (n=l) read: 

3 2 17 „ 23 2 41 2 

p u = —eu + -u — 6uv — Y2 ~2~ u v — 2~ UV ' 

(3 V = —ev + uv H — 4t> 2 — j^u 2 v + —uv 2 — tt 1 * 3 - (36) 

.1. _ Zi Zi 

We shall not present the expressions for the other RG functions here, as we are going 
to study only the fixed point equations. 
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Looking for the solutions of the fixed point equations for functions (p5|), (|36|) one 
can show that in the one-loop approximation, in addition to the Gaussian fixed point 
u* = v* = 0, there exist two more solutions u* = 2/3, v * = and u* — 0, v* — —1/5 
and the solution u* ^ 0, v* ^ is absent 96 . The fixed point with v * < is beyond the 
region of parameters describing the diluted magnet 97 and the pure model fixed point 
u* ^ 0, v * = appears to be unstable with respect to the f-coupling (we propose to 
the reader to check this by looking at the stability matrix By(u,v) (0) eigenvalues 
at the fixed points). The corresponding plot of the functions /3 U , (3 V in the one-loop 
approximation is shown in the figure || Passing to the two-loop approximation makes 
the result even worse: only the Gaussian fixed point is present (see Fig. ^). Returning 
back to the initial problem statement one should conclude that the obtained picture 
corresponds to the absence of a second order phase transition in a d = 3 Ising model 
with a weak dilution as well as without a dilution (the absence of the fixed point 
u* 7^ 0, v* =0). Of course, this is contrary to the real situation. Let us further note 
that the behaviour obtained for the /3-functions of the model (^) in the one- and two- 
loop approximations (Figs. ^ f|) resembles those for the superconductor case in the 
corresponding approximations (Figs. [1], 0). 




Figure 4. /3-functions of the diluted Ising model f3 u (u,v), f3 v (u,v) in two-loop approximation for 
d = 3. Only the Gaussian fixed point u* = v* = survives. 

Applying, however, the resummation procedure to the series (|3"5j), (|3T)| ) in the two 
loop approximation one reconstitutes fixed points {u* ^ 0, v* — 0), {u* = 0, v * ^ 0) 
and obtains a new stable fixed point u* ^ 0, v * ^ which governs a second order 
phase transition in a weakly diluted Ising model. The picture obtained appears to be 
stable with respect to successive accounts of the higher order terms in the perturbation 
theory, when the appropriate resummation technique is applied. As we have already 
claimed above, these RG results are confirmed by other different theoretical approaches 
and correspond to the experimentally observed second order phase transition in the 
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weakly diluted Ising magnet with critical exponents differing from those of the pure 
case. In Fig. |5] we show the crossing of the (3 u (u,v) and j3 v (u,v) surfaces for the 
resummed function. The calculations have been performed by means of the Pade-Borel 
resummation technique for the resolvent series (0) of two-loop functions (p5|), fl36|) as 
described above 98 . The Gaussian (u* = v * = 0) and pure (u* = 1.3146, v* = 0) fixed 
points can be seen at the rear of the cube. The cross-section of u and v planes in the 
picture are chosen to pass through the stable fixed point {u* = 1.6330, v* = 0.0835) 
corresponding to a new critical behaviour. 




Figure 5. /3-functions of the diluted Ising model (3 u (u,v), f3 v (u,v) in two-loop approximation for 
d = 3 obtained by applying Pade-Borel resummation technique. Resummation restores the presence 
of the fixed point u* ^ 0, v* = and results in the appearance of a new stable fixed point 
u* ^ 0, v* ^ 0. 

The example we have considered above is a typical situation happening in a d = 
3 RG theory: when considered without an appropriate resummation technique, the 
RG analysis might not only give imprecise values for critical exponents but also a 
qualitatively wrong answer about the absence of a stable fixed point for a certain model, 
resulting in the absence of a second order phase transition. Now with this information 
in hand let us pass to the analysis of a model of a superconductor described in the 
two- loop approximation by the RG functions (pop, (pop. 



5 FIXED POINTS AND FLOWS IN THREE DIMENSIONS 



We shall continue here by considering the flow equations (|10|) directly for d = 3. 
We shall look for the solutions of the fixed point equations at d = 3 paying attention to 



the possible asymptotic nature of the corresponding series (0), 0l9|) . Consider first the 



equation for the uncharged fixed point U. Substituting /* = into (0) one obtains 
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the following expression for the function 0% = /3 U (u, f* = 0): 



lf , n + 8 9 3n + 14 , 
V — — u ■ 



6 12 

Solving this polynomial for the fixed point one obtains for the non-trivial u* > 0: 

n + 8 Vn 2 - 20n - 104 



(37) 



u 



*u 



3n + 14 



3n + 14 



and immediately the "condition of the existence of a non-trivial solution u* Un qualita- 
tively very similar to this, appearing in the frame of the ^-expansion technique, follows 
(see 15 ' 30 and formula ( pl|) of the present article as well): the solution exists only for 
certain values of n > n c = 24.3 ! 

0.500 q /3^ 



0.300 



0.100 



-0.100 



-0.300 



-0.500 



Pade-Borel resummation 
[1/1] Pade-approximant 
Non-resummed function. 




Figure 6. /3„-function of the uncharged model /?„ at d — 3, n = 2. 



From Fig. ^| one can see that the function 0^ (37) does not intersect the w-axis 
for any non-zero value of u for n — 2. In the 0(n) -symmetric 4 -theory at d — 3 
this situation is well known (see e.g. 94 ' "). The /^-function calculated directly at d = 3 
does not possess a stable zero for the realistic values of n, nevertheless with in the 
three-loop order the presence of the stable fixed point is restored. To avoid this artifact 
appearing in the two-loop calculation one can either resum the series for the /3-function 
or construct the appropriate Pade approximant 100 in order to perform the analytical 



continuation of (|37_D out of the domain of convergence (which is equal to zero for the 



series in the right-hand side of (|57)). Let us try both methods. Representing ( p?]) in 
the form of a [1/1] Pade approximant: 

6 = v~- 

l + B u u 



fffo" = u 1 + ^ uU (39) 



one obtains: 

n 2 + 7n + 22 3n + 14 

6(n + 8) 2 [n + 8 



18 



Table 1. The fixed point U coordinate function of n. u ' 

obtained on the basis of the [1/1] Pade approximant; u* u,Res : obtained 
by the Pade-Borel resummation. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


r. ,*U,Pade 

Lb 


1.800 


1.500 


1.269 


1.091 


0.951 


0.840 


0.750 


0.676 


u *U,Res 


1.315 


1.142 


1.002 


0.888 


0.794 


0.717 


0.652 


0.597 



and, solving the equation for the fixed point: 

pU,Pade( u *P,Pade^ _ g ^ 



one obtains: 



*u,Pade _ 6(n + 8) 



n 2 + 7n + 22 y J 

So we have obtained a qualitatively different situation. The behavior of the function 
f3^' Pade {u) for n = 2 is shown in Fig. |6| by the dashed curve. If one is interested in more 
accurate values of u* some resummation has to be applied. Choosing the Pade-Borel 
resummation technique 101 and following the scheme ( |26"D - (|29| ) one obtains for the 
resummed function /3^' Res98 : 

= u[2(l - A u /B u )(l - E{^-)) - 1], (43) 

uB u 



the coefficients A u , B u are given by (|40|), E(x) = xe x Ei(x) , where the function 

oo 

E^x) = e~ x J dte~\x + ty l 
o 

is connected with the exponential integral by the relation 102 : 

Ei(x ± iO) = -Ei(-x) =F in. 

The behavior of the function f3^ ,Res (u) is shown in Fig. |6] by the solid curve and 
the fixed point coordinate u* u,Res is obtained solving the non-linear equation: 

(3^' Res {u* u ' Res ) = 0. (44) 

The coordinates of the fixed point u* u obtained on the basis of the Pade approximation 
and the Pade-Borel resummation ( u * c/ ) Pa<fe j u *u,Res^ f QT different n are given in Table [l|. 

From this analysis we conclude that: in the d = 3 theory the Pade approximants 
(as an analytical continuation of /^-functions) may change the picture qualitatively 
and lead to values of fixed points comparable to those obtained by the Pade-Borel 
resummation technique. 

Consider now the equation for the charged fixed point C applying the above con- 
siderations to (3f, for which the expression at d = 3 reads (|I~8|): 

Tl 

Pf = -f + ^f 2 + nf 3 . (45) 
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Figure 7. /^/-function at d — 3, n = 2. 



The behavior of f3f as a function of / is shown in Fig. |7| by asterisks. Note, however, 
that in this case the function f3f even without any resummation possesses a non-trivial 
zero f* M (its value f* c ^ Dtr is given in the 2nd row of Table [|). Representing ( f45f ) in 
the form of the [1/1] Pade approximant: 

1 + 5// 

one has for A/, 5/: 

Af = ^, B f =-6, (47) 
and, solving the equation for the fixed point coordinate f* c > Pade \ 

pPade^*C,Padej = q ( 4 g) 

one obtains: 

j*C,Pade _ 6 

n + 36 

The function f3f (f) is shown in Fig. [7|by the dashed line, the coordinate f* c ^ Pade 
is given in the 3rd row of Table 0. But now the series (|45|) is not alternating and this 
results in the presence of a pole (at / = 1/6) in the approximant (£K]). Therefore, 
( f4"6|) correctly represents the function /?/(/) only for / < 1/6. Let us note however 
that for all positive n a fixed point exists and its coordinate f* M ' Pade lies within the 
limits < f* c > Pade < 1/6, where there are no poles in (f46"|). Comparing this result 
with that obtained for the uncharged fixed point one can note that the representation 
of /3f in the form of the Pade approximant does not qualitatively change the picture 
(a solution for /?/(/) = exists at d = 3 even without an analytical continuation) but 
results in a decrease of the fixed point coordinate. Contrary to the e-expansion values 
(j22"D there are no borderline values of n for the positivity of f* c . Unfortunately, we 
cannot check this result by means of the Pade-Borel resummation technique. The above 
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Table 2. The fixed point C coordinate / as a function of n. f* c > Dir ; 
obtained by a direct solution of the equation for the fixed point; f* c ^ Pade - m 
obtained on the basis of the [1/1] - Pade approximant; f* c ' e ; the e-expansion 
result with the linear accuracy in e; f ,£ : the e-expansion result with the 
square accuracy in e. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


j*C,Dir 


0.920 


0.629 


0.500 


0.424 


0.372 


0.333 


0.304 


0.280 


j*C,Pade 


0.162 


0.158 


0.154 


0.150 


0.146 


0.143 


0.140 


0.136 


r C,e 


6.000 


3.000 


2.000 


1.500 


1.200 


1.000 


0.857 


0.750 




-210.000 


-51.000 


-22.000 


-12.000 


-7.440 


-5.000 


-3.551 


-2.625 



mentioned presence of a pole in the denominator of the Pade approximant makes the 
corresponding integral representation problematic 103 . In order to find the w-coordinate 
of the fixed point C, u* c , we have to deal with a function of two variables, j3 u (u,f), 
represented by a rather short series (0). Another problem arises due to the fact that 
the function f3 u (u, f) contains generating terms (i.e. (3 u {u = 0, /) ^ 0). In order to 
perform some kind of an analytic continuation of the function of two variables one can 
use the Chisholm approximants d32|) 91, 93 . But, the presence of generating terms makes 
this choice rather ambiguous. The most reliable way in such case is the representation 
of (3 u (u, f) in the form of a resolvent series B(u, f, r) (|34]) 92 ' 93 introducing an auxiliary 
variable r, which allows the separation of contributions from different orders of the 
perturbation theory in the coupling constants. The series for B(u, f, r) then reads: 

B(u,f,T)=i3 u (uTjr) = J2 bjri, (50) 

i>o 

with the usual notation applied for the coefficients bj. Now one considers ([50]) as a 
series in the single variable r. This series can be represented in the form of the Pade 
approximant B Pade (u, /, r) as an analytical continuation of the function B(u,f,r) for 
a general value of r. In particular at r = 1 the equality holds B(u, f, r — 1) = /3 u (u, f) 
and the approximant 

B Pad6 (u,f,T=l)=tf ad6 (u,f) 



represents the initial function (3 u {u, f). In our case the expression for B(u, f, r) reads: 

B(u, f, r) = r(6i + b 2 r + 6 3 r 2 ), (51) 



where: 



h = -u, b 2 = ^^u 2 - 6uf + 18/ 2 , 
6 

3n + U 3 2n + 10 2 71n + 174 2 , 



w + 5 u 7 + — uf 2 - (7n + 90)/ 3 . (52) 



12 3 " 12 
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Figure 8. Intersection of the function j3 pade (u, f) at d = 3, n = 2 with the plain / = f* c '. Pade [ n 
two-loop approximation. 

Table 3. The fixed point C coordinates u * c > Pade obtained on the 
basis of the [1/1] Pade approximant for the "resolvent" series as a 
function of n. CI : the unstable fixed point; C2 : the stable fixed 
point. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


CI 


0.184 


0.181 


0.179 


0.177 


0.175 


0.175 


0.176 


0.179 


C2 


3.309 


2.457 


1.781 


1.150 


0.473 


0.369 


0.305 


0.256 



Representing the expression in brackets in the right-hand side of (|5lD in the form of a 
[1/1] Pade approximant we have: 

B p «*(u,f,T)=Tb 1 )+ J %" fT , (53) 
1 + B uJ r 

where 

a h h -b 3 

Au ' f = vr^ s -' = ir (54) 

Let us note here that the function B(u, f,t) obtained in this way as the approximant 
for the function of two variables /3 u (u,f ) obeys certain projection properties in the 
single-variable case: substituting / = or u = into Q53D one obtains the [1/1] Pade 
approximant for /3%(u) or the [0/1] Pade approximant for f3 u (u = 0,f). Finally the 
expression for f3 u (u, f) approximated in such way reads: 

ATW)= ^i^-tH (55) 

Substituting into the equation for the fixed point /3 u (u* c ,f* c ) = the value for the 
coordinate f* c = f* c > Pade (OT) one obtains the non-linear equation for u* c ' Pade : 
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(3^( U j = r c ' pad4 ) = o. 



(56) 




Figure 9. /^-functions of the model of superconductor (3 u (u, /), /3/(it, /) in two-loop approximation 
for d = 3, n = 2 obtained by Pade analysis for the resolvent series. The stable "charged" fixed point 
C2 with coordinates u* = 2.457, /* = 0.158 as well as the unstable fixed point CI 
u* = 0.181, /* = 0.158 are seen on the front side of the cube. Gaussian (u* = /* = 0) and 
"uncharged" (it* = 1.500, /* = 0) fixed points are located on the rear of the cube. 

Solving ( |56"D with respect to u one obtains the values u * c > Pade given in Table [| The 
intersection of the function (3 Pade (u, f) (|55|) with the plane / = f* c > Fade is shown for 
n = 2 in Fig. The first fixed point (CI) given in the 2nd row of Table |3] turns out to 
be unstable, while the fixed point C2 is stable in the case n = 2 we are predominantly 
interested in. 

The resulting picture of /3-functions surfaces is shown in the Fig. |9] . The Gaus- 
sian and uncharged fixed points may be seen at the rear of the picture, whereas the 
intersection of the u- and /-planes has been chosen in the picture to cross the stable 
fixed point C2. The unstable fixed point CI is seen as well. 

The crossover to the asymptotic critical behavior is described by the solutions of the 
flow equations ( |10D with the initial values of u(£ ) and f(£o) at £ = £q 10A . Substituting 
the analytical continuation of the /3-functions in the right-hand side of (|T0D, with the 
Pade approximants (|46|) , (|55|) we get the following system of differential equations: 



L dl ~ J 1 + B f f ' 1 j 

du _ 1 + Auj . . 

dl ~ 1 + B U / 

where Af, Bf and A u j, B u j are given by ([47]) and (|5~4"1), respectively 
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Figure 10. Flow lines for the case n = 2, d = 3 given by equations (57), (|58|). Fixed points G, U, 
CI are unstable, fixed point C2 (shown by a box) is a stable one (for further description see text). 



Solving equations (57Q, (|58|) numerically one obtains the flow diagram shown in 
Fig. [10] for the case n = 2. The space of the couplings is divided into several parts by 



separatrices (thick lines in Fig. 10) connecting the fixed points. Besides the Gaussian 
(G) there exist three fixed points, one corresponding to the uncharged (U) and two 
others corresponding to the charged (CI, C2) cases. The fixed points G, CI and U are 
unstable (solid circles in Fig. [IT]) and the fixed point C2 is stable (shown as a solid box 
in Fig. UCj). Several different flow lines are shown in Fig. |H| They can be compared 
with the corresponding flow picture obtained by a direct solution of the flow equations 
for the two-loop /3-functions expressed by the third-order polynomials in couplings u, 
f (|18D , (|19D (see Fig. 2a in 30 ). There, one can see that no stable fixed point exist and 
furthermore that the fixed point U is absent. Comparing Fig. |lCj and Fig. 2b from 30 
one can see how an analytical continuation of the /3-functions (|10D , (|18|), done only 
partly in 30 and performed here in the form of Pade approximants restores the presence 
of the fixed point U (unstable) and leads to the appearance of a new stable fixed point 
C2 for the charged model. The coordinates of the fixed points U, CI, C2 are given in 
the corresponding rows of Tables [I], |2], |3] and for n = 2 they are equal to: 



U : u* = 1.500, /* = 0, 

CI : u* = 0.181,/* = 0.158, 

C2 : u* = 2.457, /* = 0.158. 



6 CRITICAL EXPONENTS 



Values of the critical exponents can be determined by the fixed point values of the 
^-functions defined on the basis of renormalizing Z-factors (|14]) - fllTD by: 



d = fidliaZi/dfi, 



(59) 
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where the derivative is taken at fixed unrenormalized couplings. The expressions for the 
(■-functions related to the order parameter and the temperature field renormalization 
in the two- loop approximation read 30 : 

C^-3/ + ^ + (H^, (60) 

s 6 12 

2(n + 2) (5n + l) 2 



3 



uf- K ' ' f, (61) 



U = ^f + nf 2 . (62) 
o 

If there exists a stable fixed point, the critical exponent of the correlation length u, of 
the order parameter susceptibility 7 and of the specific heat a will be given by: 

u = (2-Q)-\ (63) 

7 = (2- 0^(2 -Q, (64) 

a = (2-C)- 1 (e-2C), (65) 

V = Q, (66) 

where ( u = (v> — ( t . The exponents (|63|) - (|6~6f) are related by the familiar scaling 
laws. From the analysis above it follows that the charged fixed point C2 is stable and 
this results in values for exponents fl63|) — (]65f) that are different from the values of the 
uncharged fixed point U, i.e. they are not given by the A He values as it is sometimes 
stated (see e.g. 16, 59, 25 ). 

Recently, an interesting consequence of the existence of a stable charged fixed point 
(C2) has been observed 62 . According to the charge renormalization ([|) the /^/-function 
reads 

P f = f (e-UM). (67) 

Thus at a fixed point with a nonzero /* the value of the gauge field (-function is given 
by Q = e exactly. This means that the penetration depth A and the correlation 
length £ are proportional and the temperature dependence follows a power law with 
the exponent u 62 . This is not the case, though, at the fixed point with /* = 0. There we 
have £4 = (each loop contribution to the (^-function contains at least one /-factor) 
and the penetration depth behaves as A ~ £( 2 - £ )/ 2 _ Thus one would have two different 
critical length scales. 

Now, trying to obtain numerical values of the critical exponents on the basis of the 
values of fixed point C2 coordinates f* c ' Pade ^ u *C2,Pade gj ven j n Tables 0, 0, in order 
to be self-consistent let us perform the same type of an analytical continuation for the 
(-functions series as those which have applied to the /3-functions (|T8|), (|i"9D. So, after 
introducing the auxiliary variable r let us represent functions (^) - QB5]) in the form 
of resolvent series in r and then chose the [1/1] Pade approximants for these series, 
which at r = 1 will give us the analytical continuation of the series requested. Thus, 
the expression for the critical exponent <f>(<j)= \y, 7, a)} reads: 



1 + B d 



The expressions for the coefficients A^, in ( j68| ) read: 

A^af + B^- B^ = -af/4\ (69) 
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and a$ are to be determined from the resolvent series in r: 

= E^M.=i- (70) 

i>0 

Substituting (|60|) and (|6l|) into (|63|) - (|65|) and representing (^) - (^) in the form of 
(|70|) one finds: 

4 0) = 1/2, 

a« = (n + 2)/12 M -3/2 /, 

= (n 2 -n- 6)/144 u 2 + (71n + 138)/48 / 2 + (n + 2)/12 it/, (71) 
af = 1, 

a« = (n + 2)/12 it, 

a^ 2 ) = (n 2 - 2n - 8)/144 m 2 + (5n + l)/4 / 2 + 5(n + 2)/24 uf, (72) 
a (0) = 1 

ai 1} = -3(n + 2)/12 u + 9/2 /, 

= (-3n 2 + 3n + 18)/144 m 2 - (71n + 138)/16 f - 
(n + 2)/4uf. (73) 

Now considering the case n = 2 and substituting coordinates of the fixed point C2 
(f*c,Pade = 15g ^ M *C2,Pade = 2.457 ( see Tables g g) into ([7Tp - ([73]) one obtains for 

the critical exponents (|63|) - (|66|) 105 : 



1/ = 0.86, 7 = 1.88, (74) 
a = -1.14, 77 = -0.19. 

The application of the Pade approximants for the analytical continuation of the 
functions may result in the appearance of poles in these functions. If the pole is located 
in a region of the expansion parameters which has no physical meaning, e.g. a negative 
coupling u or /, the analysis is not complicated. This is the case for the /3-functions 
in the region of couplings less than the fixed point values. For the ^-functions, how- 
ever, considering the non-asymptotic behavior (and thus being far from the stable fixed 
point) one passes through a region of couplings where the Pade approximation for the 
^-functions becomes ambiguous which results in the appearance of a pole. Therefore, 
studying the crossover behavior in the next subsection we will still keep the polyno- 
mial representation for (^-functions instead of the Pade approximants. Then for the 
asymptotic values of critical exponents one gets: 



v = 0.77, 7 = 1.62, (75) 
a = -0.31, 7] = -0.10. 

The comparison of the values (|75|) and ([76[ ) shows a numerical difference of 15% in v 
and 7 and a considerable increase of a but there is no qualitative change (e.g. the 
sign of the specific heat exponent remains the same). This should be compared with 
values given by others, v = 0.53 and 77 = — 0.70 62 and 77 = — 0.38 61 . Since for the 
conventional superconductors the experimentally accessible regime lies in the precritical 
region further away from T c , let us now discuss some non-asymptotic quantities such 
as effective exponents and amplitude ratios. 
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7 AMPLITUDE RATIO FOR THE SPECIFIC HEAT 



One of the most interesting measurable quantities is the specific heat. Near a 
second order phase transition, asymptotically it follows a power law 

C ± = — \t\~ a + const. (76) 
a 

where ± indicates the specific heat C and its non universal amplitude A above and 
below T c . The amplitude ratio A + / A~ found from the ratio C + (t + ) /C~(t~)o after 
subtracting the non singular background value constitutes a universal quantity at T c 
depending only on the dimension and the number of components of the order parameter. 

The calculation of this ratio can be extended to the non asymptotic region 106, 107 
resulting in a temperature dependent measurable quantity. This also tests the descrip- 
tion of the nonasymptotic behavior by a certain flow in the interaction space of the 
Hamiltonians, discussed for the effective exponents. The starting point in the calcula- 
tion is the renormalization group equation for the specific heat 



d n d n d / <9 N 



C ± (t,uJ,») = t i- £ B(uJ), (77) 



where the inhomogeneity B comes from the additive renormalization. The formal 
solution reads 



C ± (t,w, f,fi) = fi e exp 



f(e- 2( n u(x) — 
J x 



exp - (e- 2( v (x)) 



dx 



x 



(7E 



The amplitude ratio is most easily calculated by choosing the same value for the flow 
parameter both above and below T c , which means t + = —2t~ 108 . We then recover the 
asymptotic expression found in 107 



A+ 
~A~ 



,Bv + F + a 
Bu + F-a 



(79) 



where the functions B and F^ 1 are taken at the fixed point. For this function we use 
the lowest order result known from the \l/ 4 theory neglecting the coupling to the gauge 
field; B = 2n, F + = —n and F~ = 12/u* — 4. Then we have for n = 2 



A+ „ 2v-at 

— = 2 a — . 80 

A- 2v -2a + Qa/u* v 1 

In the Table ^ we have collected the values obtained for the different fixed points. 
It is interesting to note the reasonable estimate for this ratio presented in 109 for the 
superfluid phase transition. The authors have found a value of A + j A~ = 1.067, which 
is surprisingly near the value obtained in the stable charged fixed point using both 
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Table 4. Asymptotic values for the specific heat 
amplitude ratio at various fixed points. The exponents 
in the first two lines correspond to the procedure 
leading to (|75f) . The third and the fourth lines 
correspond to the procedure leading to ([761). 



F P 


A + /A- 


ut 


7 / 

£/ 


ti* 


U 


1.81 


-0.33 


0.72 


1.500 


C2 


1.07 


-1.14 


0.86 


2.457 


U 


0.78 


0.15 


0.62 


1.500 


C2 


1.06 


-0.31 


0.77 


2.457 



calculation schemes, although the exponents are very different. We have not, however, 
taken into account changes in the scaling functions due to the coupling /. 

In comparison with experimental results 18 , the amplitude ratio of the Gaussian 
n- vector model without a coupling to the gauge field A + / A~ = n/2 3 / 2 has been used 
since the dimension of the order parameter is unclear. Later on this expression for the 
amplitude ratio has been calculated for cases other than the isotropic symmetry. This 
leads to a dependence of the ratio on the higher order couplings 110 . 

8 EFFECTIVE EXPONENTS 

Effective exponents are usually defined by the logarithmic temperature derivatives 
of the corresponding correlation functions (see e.g. 26 ' m ). They can be found from 
the solutions of the renormalization group equation for the renormalized vertex func- 
tions. The effective exponents contain two contributions, one from the corresponding 
^-functions now taken at values of u(£), f(£) of the particular flow curve considered 
("the exponent part"), and one from the change of the corresponding scaling function 
("the amplitude part"). For the analysis below, we neglect the latter contributions 
since we expect them to be smaller than the differences for the fixed point values of the 
exponents coming from the different treatments discussed before. Thus we have: 

v=(2-Cv(e))-\ (81) 
7 = (2-C,W)- 1 (2-gW), (82) 
a = (2-Ue))-\e-2( u (e)). (83) 

The flow parameter £ can be related to the relative temperature distance T c by the 
matching condition t(£) = (^o" 1 ^) 2 , with £ the amplitude of the correlation length. 

We have computed these effective exponents, see Fig. |11| - Fig. [13], along the flow 
lines shown in Fig. |H] by inserting 112 values of the couplings u(£) and f(£) into Eqs. 
fl8"T|) - (B3). For the separatrix 1 we have started with initial conditions leading to a flow, 
which does not coincide with the fixed point CI but slightly misses it although the flow 
curve does not differ from the separatrix within the thickness of the lines shown in Fig. 
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Figure 13. Effective exponent a for the flows shown in Fig.[10] (for further description see text). 

p~0| . For the curve number 4, we start somewhat further away from the Gaussian fixed 
point G leading to the initial values of the effective exponents between their Gaussian 
values and their values for the uncharged fixed point U. Note that when the values of 
the effective exponent 7 for the uncharged fixed point U and the charged fixed point 
CI are the same within the accuracy given by the scale of the figure. 

Note that when the coupling / to the gauge field fluctuations is small (i.e. for 
extreme type-II superconductors) the RG flow passes very closely to the uncharged 
fixed point U and the effective exponents, and within some region of temperatures, 
they coincide with those of the uncharged superfluid liquid. In this region the effective 
Hamiltonian (|2|) may be considered as that of a superconductor in a constant magnetic 
field neglecting magnetic field fluctuations. Recently for such a model it has been 
shown that near the zero-field critical point the singular part of the free energy scales 
as F S i n g ~ \t\ 2 ~ a J r (B\t\~ 2u ) with v being the coherence length exponent 113 . 



9 CONCLUSIONS 

Does the above account give a definite conclusion about the order of the phase 
transition occurring in a model of the superconductor minimally coupled to the gauge 
field? First of all one should keep in mind that such an answer may be given in an 
analytical theory framework only by obtaining an exact result or a rigorous proof. Here, 
the problem has been treated by a perturbation theory approach and the account of 
the influence of fluctuations on the order of the phase transition is studied within the 
field theoretical RG technique. We have shown that remaining within this approach 
one may get an answer about the second order phase transition occurring in the above 
mentioned model. 

The main point discussed in this context is whether the equations for /3-functions 
possess a stable fixed point or not. The absence of the stable fixed point is often 
interpreted as a change of the order of the phase transition (caused by the presence of 
the magnetic field fluctuations) and the evidence of the fluctuation-induced first-order 
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phase transition. This change of the order of the phase transition (being of the second- 
order in the absence of the coupling to the gauge), however, is confirmed only by the 
perturbation theory calculations in low orders 15 , (see 30 and the references therein as 
well). 

Applying a simple Pade analysis to the series under discussion 114 we have shown 
how one can recover a stable fixed point in the RG equations. In the case of one cou- 
pling, such an approach gives a qualitatively correct picture of the phase transition and 
restores the presence of a stable fixed point ( 94 , see formulas (|38|) , ( f4~2"D of this article 
as well). The same occurs in the case of two couplings: for n = 2 the "uncharged" 
fixed point U (having coordinates f* u > Pade = .158, u* u,Pade = 2.457) appears to be 
stable, which leads to a new set of critical exponents. We should note, however, that 
the pair correlation function critical exponent rj, calculated by familiar scaling rela- 
tions on the basis of sets of values ( [75] ) or ([?J]), remains negative which agrees with the 
result of 61, 62 ' 54 . Being calculated only in a two-loop approximation with the applica- 
tion of the Pade analysis, these values for the critical exponents are to be considered 
as preliminary ones. The main point we wish to make is that within the framework 
of the renormalization group analysis for the superconductor model there remains the 
possibility of a second-order phase transition characterized by a set of critical expo- 
nents differing from those of 4 He. Another important task could be to calculate the 
nonasymptotic specific heat in order to compare with experiments within the region of 
the crossover to the background. 
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